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Exclusive baryon pair production through the radiative return is a unique tool to determine the 
electric and the magnetic form factors in the timelike region over a wide range of energies. The decay 
of A-baryons into p + ti~ can be used to measure their spin and to analyze spin-spin correlations 
in the AA system. We evaluate the spin dependent hadronic tensor, which describes this reaction, 
calculate the corresponding cross section for e'^e~ — » AA and e^e~ — > AA7 with subsequent decays 
of A and A and present a compact formula for the correlated distribution. We demonstrate that 
the relative phase between the electric and the magnetic form factor can be determined in this 
analysis by considering single spin asymmetries and correlated decay distributions. The reaction is 
implemented in the event generator PHOKHARA. It is demonstrated that the relative phase can 
be measured with the event sample presently accumulated at B- factories. 

PACS numbers: 13.66.Bc, 13.88.+e, 13.40.Gp, 13.30.Eg, 14.20.Jn 
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I. INTRODUCTION 

Measurements of the proton-antiproton production 
through the radiative return, as suggested in have 
lead to a determination of the proton electric and mag- 
netic form factors, Ge and Gm, in the time- like region 
over an energy which ranges from threshold to nearly 
3 GeV [2] ■ The study of the angular distributions allows 
to separate the contributions from |Ga/P and \Ge\'^ to 
the cross section. In the absence of longitudinally po- 
larized beams and without measuring of the polarization 
of the produced baryons, the relative phase between Gm 
and Ge does not enter the predictions In this paper 
we demonstrate that AA- production through the radia- 
tive return can be used to fully determine modulus and 
relative phase A(/) of A electric and magnetic form factors. 
We present a compact form for the hadronic tensor 



(AA|J-^|0) (OIJJIAA) 



J 



including its spin dependent part and analyze single-spin 
asymmetries and spin-spin correlations for AA produc- 
tion through e"'"e~ annihilation and through the radiative 
return. The self analyzing decays A —^ pir^ and A —>■ ni:^ 
can be used to directly translate the spin asymmetries 
and correlations into predictions for the decay products. 
The spin component normal to the production plane is 
a measure of sin(A0), while spin-spin correlations can 
be used to resolve the remaining twofold ambiguity and 
determine the relative phase completely. 
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In addition to this theoretical study we construct an 
event generator (an extension of the event generator 
PHOKHARA [U, 1 i, i, i, 0, i, i), which includes the 
spin correlated decays of A and A. We adopt a simple 
model for the form factor and demonstrate, that the de- 
termination of the phase is indeed feasible. Once required 
by an experiment, the approach can be trivially extended 
to the full baryon octet. 



II. 



THE METHOD 



The reactions e"'"(pi)e (7)2) A((7i, 5i)A(g2, S'2) and 
e+(pi)e-(p2) A{qi,Si)A{q2,S2h{k) (with photon 
emitted from one of the initial states) are both in Born 
approximation fully described by electromagnetic current 
with two independent form factors 



(1) J tL ^ -ie ■ u{q2, S2) 



4toa 



[7^,^] «(gi,5i), (2) 



related to the electric Ge and magnetic Gm form factors 
by 



Gm — Fi 
with T 



Ft GE = Ft + TFt 



(3) 



and Q ^ qi+ q2. 



The qi, Si, i — 1,2 are the four momenta and the spin 
four vectors of A and A with the properties St.qt = 
0, Sf — ~1. We will not consider the initial lepton polar- 
izations and denote their four momenta by pi (positron) 
and p2 (electron). The resulting squared amplitudes, 
which describe the two processes under consideration, are 
then given by 



(4) 



2 



where 0(1) stands for the process without (with) photon. symmetric (H^j^) parts defined as 
The hadronic tensor, which is identical for both cases, 

and defined as ^ ^ 

^ g A = - {JtiJl + JIJv) , H^^, = - {J^iJl ~ JIJu) ■ (6) 

J^p.L' = Jfj^Ju = ti^^ + H^^ , (5) ^ ^ 

is naturally splitted into symmetric [H^^) and anti- They are equal to 



^-^/L — i£|LL ^(5.S'2 {Qf_n >5'i^}+ + Q.Si {Q^, S'2iy}+ — >5'i.S'2 Qv — {<5'i^, 82,^}^ 
\Gm?{^ S1.S2 - Q.Si Q.S2) g^, + ( ^Zur% Q-^i ^-^2 + (\Gm\^ - ^\Ge\^) S1.S2 ) 

mA(-r-l) 



+lBMi^M£M£ll(^Q,S, {q^, - Q.5i {q,,S2A+) + 'Zfi^'^if {^P,s, Q^q-' {Sf + si), + , 



(7) 



^hA^^M£^(q,S2 {q^, S,A-~Q-Si {g^, ^2.}-) 



» ''^l^ir-if i^PfS,. Q^q^ iS( + Sl),q,A^ + z J^^Le^^^^ QPq, (Q.^^ - Q.S,) , (8) 



where g = (92 - gi)/2, {a^, 6y}± = a^foi. ± fo^a^ and 
i/^^, = 2|GAfP(gpQ.-.g^.g2) 
-;^(|GMp-i|G^p)g,,g. (9) 

is the unpolarized hadronic tensor. 

In the first step, for illustration and further reference, 
we shall discuss baryon production in electron-positron 
annihilation without photon radiation. The application 
to the radiative return will be presented subsequently. 

The lowest order leptonic tensor, when one neglects 
the electron mass, is equal to 

L% = ^(-2sg^'' + 4K,p2^}+) 



= diag(0, 1,1,0) , (10) 

s 

while the explicit form of L^^^^ (describing the radiative 
return), together with definitions and notation, can be 
found for instance in . In both cases the leptonic 

tensor is symmetric (apart from small imaginary part 
coming from the radiative corrections), hence the asym- 
metric part of the hadronic tensor does not contribute to 
the cross section 

da{e+e- ^ AA) - ^^L%H'^-^d<^>2{pi+P2;qi,q2) , (11) 

which can be written in the following compact form 



;| g.p r\GM\'-R_eiG,rG',) ^^g^ p.52 - Q.52 p.^l) + 2 S |Gm P (« + 4 p.^l P.S2) 



Q.Si Q.S2 



■Th (\Gm? \\Ge\^ (8 {q.v? + 2 . mi - 4)(5i.52 - 1) - q.v sp,,, gV(5r + Sl^)p\ 



(12) 



where p — {p2 — pi)/2 and s — (pi +^2)^ — ■ This re- suit becomes more transparent, if we use the spin vectors 



3 

of the A and A in their respective rest frame 



L%H^'^ = 4^2a2| |GmP (l + cos^ + i|G^,p sin^ + '"''^^^^ ^^29-^) [SI + 5|) 

+ {k\GE?-\GM?)sn?e-^ ^|^l-(i|G£|2sin2 0A-|GM|2(l + cos2 0A)) , (13) 

I 



where S\ {S\) is the zth component of the unit vector 
pointing into the direction of the A (A) spin, calculated 
in the A (A) rest frame, and Oj^ is the A polar angle 
in the e^e^ center of mass frame. The y-direction is 
perpendicular to the production plane with the positive 
direction defined through Cy = ee+ x e^, where ee+ and 
are the unit vectors pointing into the positron and A 
momenta directions respectively (see Fig. 1). 



Y 




FIG. 1: (color online) . Definition of the three reference frames 
used to obtain Eq. (I13|l . 



The linear polarization is a consequence of the relative 
phase between electric and magnetic form factors and 
points into the direction normal to the production plane 
[111. [T3| (and references therein). The correlation coeffi- 
cients Cx,z,Cx.xiCy^y,Cz,z SuYB lu gcucral non vanishing 
even if the relative phase between Ge and Gm is zero or 

TT. 

The absolute values of the electric and magnetic form 
factors can be measured without measuring the polar- 
ization of the baryons via the Roscnbluth method, i.e. 
by analyzing the angular distributions of the baryons. In 
previous investigations it was shown JTlj that the relative 



phase can be determined by using longitudinally polar- 
ized beams. From Eq. p^ it is obvious that a measure- 
ment of 5^, which can be performed for A, as discussed 
below, determines that phase modulo tt. The remain- 
ing twofold ambiguity can be eliminated through a study 
of the correlation between 5*^ and 5*^ (or and Sj^ ), 
which is proportional to Re{GMG*^). 

The measurement of the subsequent two body decays 
of A ^ 'K~p, TT^n and its charge conjugate allow for a 
straightforward spin analysis of the decaying As. The 
neutral modes are experimentally more demanding. In 
the remaining part of this work we will therefor refer only 
to the charged modes even if the theoretical description 
of both is identical. The decay distribution, when the 
spin of the proton is not measured, is proportional to 

i?A = 1 - "A ■ n^- , (14) 

where n^- is a unit vector in the direction of the tt^ 
momentum. The factor ua — 0.642(13) characterizes 
the analyzing power of the angular distribution. For the 
antiparticle — ~o:\. 

To combine production and decay, one sums over the 
intermediate spins of the A and A and uses the fact that 
the sums over the A polarizations can be performed eas- 

iiy, 

J2Sa = 0, Y^^ASi^S^'- (15) 

pol. pol. 

The spin vectors in the Eq. fT^l) are then effectively re- 
placed by S^A — ttA"-x- and S^a ^ -~a]^n„+. 

As a result the relation between cross sections of the 
reaction e+e~ AA and the reaction e+e~ — > A(^ 
7r+p)A(— > 7r~p), when the spins of the final proton and 
antiprotons are not measured, is simple. Using the nar- 
row width approximation one finds 

da (e+e~ — > A(— > 7r+p)A(— > n^p)) = 

da (e+e" AA) (5a^a ^ TaA"-7rT) 

X d^2{qi;Pn+,Pp)d^2{q2;p^-,Pp) 

X Br(A 7r+p)Br(A -R-p) , (16) 

where is the two body phase space normalized to 1, 
ttA is the asymmetry parameter of the A and n^+(njr-) 
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is a four vector, which in the A (A) rest frame is 
equal to (0,n^+) ({0,n^~)), with (n^-) being the 
unit vector in the direction of (Ptt-)- The symbol 
i^A A ~^ TQ^aJt-ttT ) indicates that Sa has to be replaced 
by -ttAn^- and % by aAn^+ . 

Let us now present a qualitative discussion of the ob- 
servables for the radiative return. In [1] it was shown that 
in the case of single photon emission from the initial state 
leptons, by choosing as reference frame the hadronic rest 
frame with the z-axis opposite to the photon direction 
(we will call that frame RF from now on), one gets (for 
^ s, which is fulfilled at B-factories) the following 
simple approximate expression for the leptonic tensor 



which is sufficient for a qualitative discussion. The y-axis 
is chosen in the e~^e~ cms frame as Ey = eg+ x e^, where 
eg+ and are unit vectors pointing the directions of the 
positron and photon momenta, respectively. It remains 
unchanged after the transformation from the e^e~ cms 
frame to the hadronic rest frame described above. Here 
2/1,2 = ^^57~(1 "F cos9y), 9j is the photon polar angle in 
the e+e~ center of mass frame and i,j — 1,2,3 in L'^. 
Other components of L'-' do not contribute to the matrix 
element. The close similarity with Eq. ([TU)) is obvious. 
As a result 



(47ra)2 (1 + cos2 6'..y) 

Q^ym 2 



diag(l,l,0) , 



(17) 



(l + cos2 6^, 
iQ yiy2 V ~/ 



(1+C0S2 Oa) 



^\Ge\' sin^ 9a - «A^^^^(^f^a) sin(20A) «- " <+) 



<i^\GE? + \GM\''' 



{^\Ge\'-\Gm\') sin' Oa nl,nl_ +al{^\GE?sm'eA-\GMmi + cosHA)) K^K- ^ ^ 



with the angle 6a being defined now in RF frame and 
(n^-) in the A (A) rest frame. Hence after measur- 
ing \Gm\ and \Ge\ with the method proposed in [1] and 
successfully used by BaBar for the proton form factors 
(no information on the spin is necessary in this case) , 
one can measure the combination cka sin(A0), where A(j} 
is the relative phase of the two form factors. This can be 
achieved for example by measuring the asymmetry 



± _ da{a^ > 0) - da{a^ < 0) 
~ da{a± > 0) + da{a± < 0) 



where a+( ^ = sin(20A) ''^^+{tt-)- Since Q!a is known with 
good precision [l^, sin(A(/)), and thus A(j), is determined 
up to a twofold ambiguity, which can be resolved by mea- 
suring the correlations between and n'^+- In prin- 
ciple, by studying the double spin correlations one can 
measure both the sign of cos(A0) and ua- However, the 
expected number of events at B-factories amounts to a 
few hundred to a few thousand distributed in the energy 
range 4m^ < < 10 GeV^, depending on the accumu- 
lated luminosity and the angular cuts in the experimental 
analysis (see Section IIVI) . The measurement of oa thus 
cannot be competitive to the existing measurements [isj . 



III. THE IMPLEMENTATION INTO THE 
EVENT GENERATOR PHOKHARA 

To allow for quantitative studies we have implemented 
theprocess into the event generator PHOKHARA P, H, 
0, S H 0, H, • This implementation will be available 
within the version PHOKHARA6.0. We use the leading 
order matrix element only, since the expected cross sec- 
tion is small and correspondingly the expected statistical 
precision is low. The radiative corrections are expected 
to be of the order of a few percent p?|. The generation 
of events is based on the method used in the generator 
TAUOLA l3|, performing first the generation with the 
unpolarized cross section da{e'^e'^ — > AA) and generat- 
ing in the second step the proper distributions of the de- 
cay products by evaluating the ratio of the cross sections 
with and without spin effects. Standard tests (see for ex- 
ample [6]), where the generation results were compared 
with existing analytic formulae, assure the technical ac- 
curacy of the implementation at the level of 10^"*. 

The cross section (T(e^e~ — s- AA) was measured at one 
energy only, and with poor accuracy. No form factor 
has been extracted. Thus one has to relay on theoret- 
ical assumptions. For the isospin and SU(3) structure 
of the electromagnetic form factors of the baryon octet 
we adopt the model ansatz suggested in reference p^. 
There all form factors are given in terms of co-, 4>- and 
p-dominated contributions. Indeed, the p piece does not 
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contribute in the case of AA production. The adopted 
model predicts real form factors Gm ^^nd Ge and thus 
zero relative phase A^. The size of the relative phase 
between Gm and G e does not alter the cross section and 
thus till now no experimental information is available on 
A(/). To show the effects of a possibly large relative phase 
we add an ad hoc phase in the analysis presented in the 
next section. 



IV. FEASIBILITY STUDY FOR B- FACTORIES 
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FIG. 3: (color online). The asymmetry Ay (see Eq.(IT5J for 
definition) for the maximal relative phase (Acji = tt/S) be- 
tween the Gm and Ge form factors. 



As. More refined strategies of this measurement are con- 
ceivable if the asymmetry is studied as a function of 0a- 
Furthermore, the sample is doubled by considering A and 
A decays. 

The combination Re[G mG*^) (or cos(A(/))) can be ex- 
tracted by measuring the asymmetry 



FIG. 2: (color online). The predicted differential cross section 
of the reaction e^e~ A(^ 7r~p)A(^ ■k'^p)^ 



The predicted cross sections as functions of Q^, with 
no cuts on the final state particles (open circles) and with 
angular cuts roughly corresponding to the angular range 
covered by the BaBar detector (filled circles), are shown 
in Fig. [21 The integrated cross sections, for \fQ^ values 
from threshold up to 3 GeV, are equal to 18 fb when 
no angular cuts are applied and 1.3 fb, if we restrict the 
polar angles of the produced charged particles between 
30 and 150 degrees. So typically at B-factories we expect 
about 130 events per 100 fb^^ accumulated luminosity. 
With the already accumulated luminosity at BaBar (390 
fb"^) and BELLE (710 fb^^) the A form factors can be 
studied with a decent accuracy. 

As already said, the adopted model predicts almost 
zero asymmetry . However, a priori a much larger rel- 
ative phase between the Gu and Ge could be present. In 
FigOthe predicted asymmetry is shown for the maximal 
relative phase between the form factors, ft is clear that in 
this case the measurement of the phase is feasible. Thus 
the experimental studies of the asymmetry will allow to 
obtain further constrains on the form factor models sup- 
plementing the information obtained by the \Gm\ and 
\Ge\ measurement through angular distributions of the 



_ (icr(a > 0) - rfcr(a < 0) 
^ duia > 0) + duia < 0) ' 
a = sin(20A) ■ (20) 

Knowing both sin(A(/)) and cos(A0) the relative phase 
is fixed. The predicted asymmetry Axz for A(f) — tt and 
a B- factory energy is shown in Fig. [H Measuring other 
spin correlations or respective asymmetries can in prin- 
ciple allow for the measurement of a\ and independent 
measurements of \Gm\ and \Ge\- On the one hand, due 
to the limited statistics, correlations studies alone will 
not be competitive with the existing measurement of cka 
or with the extraction of \Gm \ and \Ge \ by means of 
the method proposed in On the other hand, for each 
value of the distribution, which is differential in cos 9a 
and in the variables describing the decay distributions of 
A and A, depends on the three quantities \Gm\, \Ge\ 
and A(j) only. A more refined experimental investigation 
of this distribution, using the generator PHOKHARA, 
may well lead to a reasonably precise determination of 
these three numbers. 



V. SUMMARY 

A compact formula has been presented for the spin 
dependent baryon production through e~^e~ annihilation 
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FIG. 4; (color online). The asymmetry A^z (see Eg. pUjl for 
definition) for the maximal relative phase (Acji = tt) between 
the Gm and Ge form factors. 



and through the radiative return. It describes the de- 
pendence on a single spin as well as spin-spin correla- 
tions. Single spin asymmetries can be used to determine 
the relative phase between the electric and the magnetic 
form factor up to a twofold ambiguity, which can be 
resolved by analyzing spin-spin correlations. The con- 
struction of an extension of the Monte Carlo event gen- 
erator PHOKHARA allowed us to demonstrate that the 
spin dependent decay distribution of A (and A) can be 
used to perform this analysis under realistic experimen- 
tal conditions. Given enough data, the study could be 
extended in a straightforward way to SS, SS production 
with the subsequent decays — > mr^,pTT^, S~ — > mr~ , 
S*' Att" and Att^ as typical examples. 
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